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Abstract 

We report on the next-to-leading order(NLO) QCD computation of top-quark pair pro- 
duction in association with a photon at the Fermilab Tevatron RUN II and CERN Large 
Hadron Collider. We describe the impact of the complete NLO QCD radiative corrections 
to this process, and provide the predictions of the leading order(LO) and NLO integrated 
cross sections, distributions of the transverse momenta of the top quark and photon for the 
LHC and Tevatron, and the LO and NLO forward-backward top-quark charge asymmetries 
for the Tevatron. We investigate the dependence of the LO and NLO cross sections on the 
renormalization/factorization scale, and find the scale dependence of the LO cross section is 
obviously improved by the NLO QCD corrections. The K-factor of the NLO QCD correction 
is 0.977(1.524) for the Tevatron(LHC). 



PACS: 14.65.Ha, 14.70.Bh, 12.38.Bx 
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I. Introduction 



The top-quark was discovered by the CDF and DO collaborations at Fermilab Tevatron in 1995 [1, 
2J. It opens up a new research field of top physics, and confirms again the three-generation 
structure of the standard model(SM)[3], H]. Among all the elementary particles discovered up to 
now[5j[6], the top-quark mass term breaks the electroweak (EW) gauge symmetry maximally due 
to its huge mass, and the detailed physics of the top-quark may be significantly different from the 
predictions provided by the SM. But until now our knowledge about top quark's properties has 
been still limited [7]. For example, the couplings of the top quark to a photon and a Z° boson 
have not yet been directly measured [H |9], while the precise measurement of the production and 
decay of top quark may be significant in searching for new physics beyond the SM. 

In recent years there have been many works devoted to the study of the top-quark couplings. 
The studies for probing the top-quark couplings ti'y and tiZ° at hadron colliders at LO were 
carried out in Ref . |10j . the calculations for the process e + e~ — > tiZ° at LO, QCD and EW NLO 
are provided in Refs.[llj. and the one-loop SM QCD and the supersymmetric QCD effects in the 
process of 77 — > tiZ° at the ILC was investigated in Ref. [12]. The SM couplings of tiV{V = 7, Z°) 
may be modified by the new interactions and that would lead to abundant phenomena of new 
physics. For example, if the top quark was a composite object, there would be an anomalously 
large tij event rate at colliders, due to deexcitation of high-energetic top state [13]. And if there 
exists nonstandard CP violation, in particular Higgs sector CP violation, a sizable top-quark 
(weak) electric dipole moment could be induced[14j. Other relevant references[9] [151 HS1 E] 
indicate that the vector and axial form factors in the coupling of the top quark and neutral gauge 
boson V(= 7, Z°) should be probed precisely in order to find the signatures of a certain model of 
dynamical EW breaking. 

At a linear collider, it is not easy to obtain the information about the individual EW neutral 
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coupling tiV (V = Z°, 7) from the precise measurement of the top-pair production at a linear 
collider because of the hardness in distinguishing the contributions from the tiZ° and ti'j couplings. 
At a hadron collider, it is impossible to measure the EW neutral couplings via qq — > 7* / 'Z* — > it 
due to the strong interaction process qq(gg) — > g* — > it. Instead, they can be measured in QCD 
tiZ /'-/ production and radiative top-quark decays in ti events (ti — > r yW + W~bb). Each of the 
processes is sensitive to the EW coupling between the top quark and the emitted Z°-boson(or 
photon). In the work of [TS] it is concluded that it will be possible to probe the ££7 coupling at 
a few percent level at the LHC Since the LO predictions in the QCD expansion for the channels 
pp(pp) — > tiZ°( , y)+X at hadron colliders contain significant theoretical uncertainty, it is important 
to improve the theoretical prediction in order to accommodate the experimental measurement of 
the top-quark couplings. Recently, the NLO QCD correction to tiZ production at the LHC has 
been calculated in Ref . [19] . 

Our study in this work corresponds to the investigation on the production of the top-quark 
pair associated with a photon at the Fermilab Tevatron Run II and the CERN LHC in both LO 
and NLO QCD approximations. It is arranged as follows: In Sec. II we provide descriptions of the 
analytical calculations. In Sec. Ill we present some numerical results and discussions, and finally 
a short summary is given. 

II. Description of the calculation 

In the calculations at the LO and NLO of the a s expansion, we use the 't Hooft-Feynman gauge, 
employ FeynArts3.4 package [20] to generate Feynman diagrams and their corresponding ampli- 
tudes. The LO amplitudes are precessed by adopting FormCalc5.4 programs [2 lj. In the calculation 
for virtual corrections, the one-loop amplitudes involving UV and IR singularities are handled an- 
alytically by using our modified FormCalc programs, and are output in Fortran code with the UV 
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and IR "ex N-point integrals" terms remained unprocessed. The output is further processed nu- 
merically by using our developed Fortran subroutines for calculating N-point integrals to extract 
the remaining finite e- terms. In these Fortran codes the IR singularities are separated from the 
IR-finite remainder by adopting the expressions for the IR singularity in N-point integrals(iV > 3) 
in terms of 3-point integrals [22j. 

II.. 1 Born approximation 

We consider five partonic processes which contribute to the process of top-pair production asso- 
ciated with a photon at LO for hadron colliders. They are gg —>■ tij and qq — > tt^f (q = u, d, c, s) 
production channels. We take the constraint for the transverse momentum for radiated photon as 
p^P > Px\uv e -S-' Print = 20 GeV. We express these partonic reactions as 

?(pi) +q(P2) -> t(p 3 ) +t(Pi) +7(P5), (q = u,d,s,c), (2.1) 

and 

g{pi) + g{p2) -> t{ P3 ) + f(p 4 ) + 7(ps), (2.2) 

where we denote the external four-momenta by pi(i = 1, ...,5) for the partonic processes qq — > 
tt'y and gg — > ti'j , separately. There are 4 LO Feynman diagrams for partonic process qq — > 
^7 (shown in FigJTJ), and 8 tree-level diagrams for the gg — > tt7 partonic process(shown in 
Figf5]). Despite being massless for photon and light-quarks(g = u,d,s), the cross sections for the 
above partonic processes are still "infrared safe" due to our constraint for the photon transverse 
momentum. 

The expression of LO cross section for the partonic processes qq — > ti r y and gg — > ttj have the 
forms respectively as 

i i ('9-^4 color i i ^o—^ co ' or 

^is^Ei^W' *:» = isirSK« (2.3) 

spin spin 

4 
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where the factors |, | and factors |, ^ in Eqs. (12.31) come from the averaging over the spins and 



colors of the initial partons, respectively, s is the partonic center-of-mass energy squared, M. LO 
and .M^q are the amplitudes of all the tree-level diagrams for the partonic processes qq — > ££7 and 
gg — > tij respectively. In above two equations the summations are taken over the spins and colors 
of all the relevant particles in the qq — > tij and gg —>■ ££7 partonic processes. The phase-space 
elements dflf and dfl^ 9 in Eqs. (12.31) is expressed as 

dsif' = L + n - 1 ») It (4§ev • (2 ' 4) 

\ i — 3 / j — 3 

According to the factorization theorem for hard scattering processes in QCD, the LO differ- 
ential cross section for the process pp(pp) — > ££7 + X at the Tevatron(LHC) can be obtained 
by performing the following integration of the differential cross section for the partonic processes 
qq — > ££7 and gg — > ££7 over the initial partonic luminosities [see Eq. (l2.5j) ] . 

ss,cc,gg i i 



dcrio — / dx\ I dx 2 



1 + 5ij 



ij=uu,dd 

da% do\. 



X 



G i/Pl (xi,fif) Gj/p 2 (x 2 ,fi f ) -jr- + Gj/ Pl (xi,fi f ) Gi/p 2 (x 2 ,Hf) -jr- dy t , (2.5) 



where yt and yt are the rapidities of the top-quark in the proton-(anti)proton and partonic center- 
of-mass systems, respectively (y t = |ln ^ §t~f ^ ; Vt = |ln ( f«^t ) )• The direction of the z-axis 
of the hadronic center-of-mass system is defined as the orientation of incoming hardron Pi (for 
the parent process P\P 2 —* ££~7 + X), while the z-axis of the partonic center-of-mass system is 
set as the orientation of radiated parton i(or j) from Piffor the partonic process ij(or ji) — > 

da . dcr .- 

££7]. The differential cross sections and are expressed in their own partonic center-of- 
mass frames, respectively. Under a boost in the z-direction to a frame with velocity /3, we have 
Vt — Vt — tanh -1 (3 and dy t = dy t . G^y^x, /J,f)(i = u,d,s,c,g, j = u,d,s,c,g) are the parton 
distribution functions(PDFs) of (anti)proton A(= P 1; P 2 ) which describe the probability to find 
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a parton with four-momentum xpA in (anti)proton A. The partonic colliding energy squared 
s = X1X2S, where s is defined as the center-of-mass energy squared of the proton-(anti)proton 
collision, ///is the factorization energy scale. In our LO calculations, we adopt the CTEQ6L1 
PDFs[23]. 

Our LO calculation shows when we take p^Lt = 20 GeV, the LO integrated cross section for 
the ti'j production is dominated by the gluon-gluon fusion partonic channel with about 66.3% at 
the LHC , while about 99.3% is contributed by the q — q(q = u, d) annihilation partonic channels 
at the Tevatron RUN II. 

II.. 2 NLO QCD corrections 

The NLO QCD corrections to the pp(pp) — > ttj + X process are contributed distinctly by the 
following four parts: 

1. the real gluon emission partonic processes qq, gg — > ti'jg, (q = u, d, s, c). 

2. the real light- (anti) quark emission partonic processes q(q)g — > ti^q{q), (q = u, d, s). 

3. the virtual corrections at the NLO to the partonic processes qq, gg ti'y, (q = u, d, s, c). 

4. the collinear counterterms of the PDF. 

In all the NLO calculations we use the dimensional regularization(DR) method in D = 4 — 
2e dimensions to isolate the UV and IR singularities. To describe the cancelations of the IR 
singularities in our calculations more clearly, we decompose the collinear counterterms of the 
PDF, 5Gi/p(x, fif) (P = p,p ; i = g,u,u,d,d,s,s), into two parts: the collinear gluon emission 
part SG^ l p 0n \x, fif) and the collinear light-quark emission part 5G^p Tk \x, [if). Their analytical 
expressions are presented as follows. 

SG q(g)/P {x, /if) = SGf^(x, (if) + SG { ^yp\x, (if), (q = u, u, d, d, s, s), (2.6) 
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where 



sr (gluon) ( \ 

dLr q(g)/P \ x ^f) 



-f(quark) 
*q/P ' 



6G q 

6G g/P '\?>Pf) 



i(quark) . 




1 dz 

— P qq{99){ Z ) G q{g)/p( X / Z ^f)i 



Z 

1 dz 



z 



Pqg(z)G g/P (x/z,Hf), 



d,d,s,s 



q=u,u 



1 dz 



P gq (z)G q/P (x/z,ii f ). (2.7) 



The virtual corrections to the processes pp(pp) — > qq, gg —>■ ti'j + X contain both soft and 
collinear IR singularities. These singularities can be canceled exactly by adding the contributions 
of the real gluon emission processes gg, gg — > ti'jg and the collinear gluon emission part of the 
PDF counterterms 5G^ l g U y p \ The real light-quark emission processes g(g)g — > ti'yq(q) contain only 
the collinear IR singularities. It can be canceled by the contributions of the collinear light-quark 
emission part of the PDF counterterms $G q q ? g yp exactly. All of these cancelations are verified 
numerically in our numerical calculations. The explicit expressions for the splitting functions 
Pij(z), (ij = qq, qg, gq, gg) can be found in Ref.j 



A. Real gluon emission corrections 



We denote the partonic processes with real gluon emissions as 
q(pi)+q(p 2 ) ^t(p 3 )+i(p A )+^(p 5 )+g(p 6 ), g(pi) + g(p 2 ) -> t(p 3 ) +t(p 4 ) +7(Ps) + 0(pe)- ( 2 -8) 

The real gluon emission partonic process gg — > tijg includes 24 LO graphs shown in Figl31 and 
the gg — > ti'jg subprocess involves 50 LO graphs(shown in FigJl]). The figures (1)-(10) in FigJUare 
s-channel diagrams, the FigsJ3Jll)-(30) are t-channel diagrams. The u-channel diagrams for the 
gg — > ti'jg subprocess are not drawn in FigSJ but can be obtained by exchanging incoming gluons 
in each t-channel diagram in FigJH The process cc — > ti'jg contains only the soft IR singularity, 
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while qq — > tt^g, (q = u, d, s) and gg — > tt^g contain both soft and IR singularities. The soft IR 
singularities can be extracted by adopting the two cutoff phase-space slicing(TCPSS) methods [21] 
respectively. The soft IR singularities in the partonic processes qq — > tt^g, (q = u, d, s, c) and 
gg — > ti'-fg at the LO cancel the corresponding soft IR singularities arising from the one-loop 
virtual corrections to qq — > tij, (q = u, d, s, c) and gg — > tij processes, respectively. 

We split the phase-space of the gg(qq) — > tijg, (q = u, d, s, c) partonic process into two regions, 
E 6 < 5 s Vs/2(soft gluon region) and E% > 5 s Vs/2(hard gluon region). Except for the cc — > ti'yg 
process, the hard gluon region is divided into hard collinear region(HC)(— ii 6 or — i 2 Q < 5 c s) and 
hard noncollinear (HC) region(— £ 16 or — 1 26 > 5 c s), where £y = (pi — pj) 2 and (pi — Pj) 2 for the 
qq — ► ti'yg and gg — > ti'yg respectively. Then the cross sections for the real gluon emission partonic 
processes can be expressed as 



a 



ii-'j 



Ji 



(ij -> Wyg) 



9,ij 9,ij 



9,1-3 1 9M 



(ij = uu, dd, ss, cc, gg) 
(ij = uu, dd, ss, gg) 



(2.9) 



The differential cross section for the partonic processes qq — > ttjg in the soft region is given as 



da 



9m 



"2vr 



1 7 1 

«(<7l2 + 9?,a) - -(#13 + 92i) - -(#14 + #23) 

00-3 



a s T(l — e) f 4n[i 
27rr(l - 2e) 



72" + ~ + A o. 



dCTqq 



da° 

qq I UO qqi 



where da^ are the LO differential cross sections for the partonic processes ij 
uu, dd, ss, cc). The soft integrals gij(i = 1, 2, 3, j = 2, 3, 4) are defined as 



9ij(Pi,Pj) 



(27r/i r ) 2e f d^pe \ 2(pi- Pj ) p 2 



p) 



(2.10) 
(2.11) 



27r JE 6 <s a ^I/2 E 6 l(Pi -P6)(Pj- Pe) (Pt-Pj) 2 (Pj-Pef 

The explicit expressions for the soft integrals gij(pi,Pj) relevant to our calculations for the qq, gg — > 
ti'yg partonic processes, can be found in Ref. [25]. By using Eqs. fl2.10H2.lll) and the related soft 
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Figure 3: The LO Feynman diagrams for the real gluon emission partonic process qq — > ttjg (q = 
u, d, s, c). 
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Figure 4: The LO Feynman diagrams for the real gluon emission partonic process gg — > tt^yg . 
The diagrams obtained by exchanging two initial gluon lines in (11)- (30) are not drawn. 
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integral expressions, we can express the coefficients A^ q g and Af $ in Eq. fl2.10p for the massless 
gg-fusion processes qq — > ti^g (g = u, d, s) in the forms as 

8 



A s 



Af 



3' 

8 16 



<-N 



1 



P3-P4 



3 \ 1 / 2 {s 34 ,m 2 ,mf 



I (log VI + log 

3 \ P\VTL t P2 m t 



log (0-340-43) 

-KlogV + log^ 1 

3 \ P2 m t 



And for the massive cc-fusion partonic process cc — > ££7*7, we get 



(2.12) 



4 5 



0. 
16 



?9 



Pi • P2 , / v . P3 • P4 . , 

log (0-12O-21) + — Tjr-. 2 sr- log (0-34 O-43 



3 \^ 1//2 ( s 12> "^c; m c) 



A 1/2 (s 34 ,m < 2 ,m 



7 

+ 7T 



3 VA 1 /2( Sl3 , m 2 



Pi • Pa , , x . P2 • P4 , , 

log (o-i 3 o- 31 ) + -— - — _ 2 log (0-240-42, 



A 1 / 2 (s 2 4,m 2 ,m 2 ) 



2 



H^w^ log (^41)+^ 



P2 -P3 



log (0-230-32) . (2.13) 



with cr. 



_ AV^j mg.ma) \l/2/„. ^2 m 2A 



Si,- + m, 2 — m 2 ) 2 — Asumj. 



For the — > partonic process in the soft region, we have 

'256 



da 



<).<)<) 



12tt 



D2 + I6D4 |.Mf I 2 



j -y D x + I6D3 ) 2Re(.Mf ] ■ Mf) 



dOf, 



(2.14) 



where the summation is taken over the spins and colors of initial and final states, and the bar over 
the summation represents taking average over the spins and colors of initial partons, and 



Mf = Mf + \m 9 s 9 , M 9 2 9 = M 99 - l -Mf, 



(2.15) 



Ai 99 , M. 99 and A4 99 are the amplitudes for s-, t- and w-channel diagrams of partonic process 
gg — > tt'y separately, and 



Mf 



^C 99 + C 99 + C 99 j M{ 9 + Qcf - C 99 + C 99 ^j M% 9 . (2.16) 
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The color factors are expressed as 



Cf = 5 ClC2 l, CI 9 = if clC2C X c , C 3 = d ClC2C X c 



(2.17) 



where f abc and d abc are antisymmetric and symmetric SU(3) structure constants respectively, 1 
and A c are identity and Gell-Mann matrices. Ci,C2 are the color indices of initial gluons, c is the 
color index of propagator gluon and 



Di = 9g 12 + 9g 13 + 9g 24 - #34, D 2 = 9g 12 + 9g 23 + 9g u - #34, 

D 3 = 6(5-12 -#14 -#23 +#34), ^4 = 6 (#12 -#13 - #24 + #34)- (2.18) 

The cross sections for the processes pp{pp) — > ij — > tt^g + X, (ij = uu, dd, ss, gg) in the hard 
noncollinear region, a?£j, and pp(pp) — > cc — > tijg + X in the hard region, a^ c5 , are finite and can 
be calculated by using Monte Carlo method. The differential cross section in the hard collinear 
region, da^jj, can be obtained by using 

,2\ n 



da 



HC 



1 



1 + Sij 



ol s r (1 — e) / 47T/! 



_2vrr(l -2e) V, s 



, <ixi / dx 2 
e / ./o Jo 



+ y 

+ (*<->• j 



1_5s dz (\ - z 
z 



z 

dz f 1 — z 
z V z 



Pa (z,e) G i/Pl (x 1 /z,ii f )G j /p 2 (x 2 ,fx f ) 
p jj (z, e) G i/Pl (xi, /i/) Gj/p 2 (x 2 /z, n f ) 



da. 



ij 



(2.19) 



where Gi(jyp(x, /i/) is the PDF of parton i(j), and P refers to (anti)proton. Pu(z, e) [i = q for q — q 
annihilation subprocess and i = g for # — # fusion subprocess) are the D-dimensional unregulated 
(z < 1) splitting functions related to the usual Altarelli-Parisi splitting kernel [26]. They can be 
written explicitly as 



P li (z,e)=P tl (z)+eP> l (z) 
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P qq { z) =C F ~-, P' qq {z) = -C F {l-z), 



P gg {z) = 2N 



: 1 — z ^ , 
+ + z(l - z) 



P'Jz) = 0, (2.20) 



1 — z z 
where N = 3 is the color number, Cp = 4/3. 

B. Real light- (anti) quark emission corrections 

Since the LO contributions from the real light- (anti) quark emission partonic processes q(q)g — > 
ttjq(q) are at the same a s order as previous real gluon emission partonic processes qq — > ti'jg and 
gg — > ££7(7 in perturbation theory, according to the Kinoshita-Lee-Nauenberg (KLN) theorem [2 7j. 
we should consider these subprocesses too. The LO Feynman diagrams for the partonic processes 
q(q)9 ~^ ttjq(q) (q — u i d, s) can be obtained by exchanging initial (anti)quark and final gluon in 
corresponding diagrams in FiglSJ 

In order to avoid the additional IR singularity at the LO for the partonic processes q(q)g — > 
ti^q{q) (q = u, d, s) due to the radiated photon from a massless light- (anti) quark, we take a photon 
transverse _„ cut and an angle cut between the jet and photon, e.g. p?> > 20 GeV and 
7 jet > 9™jet = 3° [in the proton-(anti)proton center-of-mass system.]^. Then these partonic pro- 
cesses contain only the initial state collinear singularities induced by strong interaction. Splitting 
the phase-space into collinear and noncollinear regions by introducing a cutoff 5 C , we can express 
the cross sections for the partonic processes qg — > ttjq and qg — > tijq as 

° R (q(q)g -> ti iq (q)) = a% )g = a% g + af {q . )g (2.21) 

The cross sections in the noncollinear region, are finite and can be evaluated in four dimen- 

sions by using Monte Carlo method. The differential cross sections in the collinear region for the 



1 From the experimental point of view, we should apply angle cut #S"| e t not only to the pp — > q(q)g —> ttjq(q)+X 
processes, but also the pp — > qq, gg — > ti^/g+X processes. The subscript 'jet' in # 7 j e t represents the light- (anti)quark 
jet or gluon jet for the real light-(anti)quark emission processes or the real gluon emission processes. 
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PP ~^ q(q)9 ~^ tt r yq(q) + X, (q — u, d, s) processes, da^-y can be written as 



da 



c 



and 



a s r (1 — e) / 47r/i 



,2\ e ' 



2vr T (1 - 2e) V s 
1 dz (\ - z" 



- ) 8 r e / dxi / dx 



Z \ Z 



P qg (z, e) G g/Pl {xi/z, n f ) G q/P2 (x 2 , Hf) 



VI 



" x dz ( 1 — z 
+ — 



+ 



+ 



z \ z 
1 dz ( 1 — z 



.''2 



xi % V z 

1 dz ( 1 — z 



P gg (z,e) G q/Pl (xi, Hf)G g/ p 2 (x 2 /z,fi f )da i 
P gq (z,e) G q/Pl (x 1 /z,/j, f )G g/ p 2 (x 2 ,/i f ) 



J-2 



z \ z 



Pgq (z, e) G g/Pl (x u iif) G q/ p 2 (x 2 / z, Hf) 



da 



99 



(2.22) 



daf = da% {q q) . 



(2.23) 



where the splitting functions P g g( g q)(z,e) can be written explicitly as 
*W*> c) = P gSlfl9 (z) + e/^Jz), P gs (z) = i[z 2 + (1 - ^) 2 ], J^(z) = -z(l - z) 



1 + (I _ 

P OT W=CV^ >- P'(z) = -C F z. 



(2.24) 



C. Virtual corrections 



There are 118 diagrams for the partonic process qq — > ti^j in the SM at NLO. It involves 
self-energy (40), vertex(32), box(14), pentagon(4) and counterterm(28) graphs. For the partonic 
process gg — ► tij there are 306 diagrams at NLO in the SM, including self-energy (32), vertex(156), 
box(66), pentagon(12) and counterterm(40) graphs. Among all these graphs the pentagon dia- 
grams are the most complicated ones. We depict them in FigjS^for partonic process qq —>■ ti^j ) 
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q 9 ' 1 g 

(1) (2) (3) " (4) 

Figure 5: The pentagon Feynman diagrams for the partonic process qq — > tij (qq = uu, dd). 

and Figj6](for partonic process gg — > ti r y ). In the NLO calculations, we use the dimensional reg- 
ularization method and adopt the modified minimal subtraction (MS) scheme to renormalize the 
strong coupling constant and the relevant masses and fields except for top-quark and gluon, where 
their masses and wave functions are renormalized by applying the on-shell scheme. The total NLO 
QCD amplitudes of partonic processes qq — > ti r y and gg — *> ti r y are UV finite after performing the 
renormalization procedure. Nevertheless, they still contain soft /collinear IR singularities. 
The virtual corrections to the subprocesses qq — > tt^y and gg — > tt^y can be expressed as 

i i (r>-\4 color 

spin 

1 I ('0 7r \4 color 

K = l^f- E 2Re [ M LoM V ga ] dnf (2.25) 

spin 

where M^o an d -Mio are the LO Feynman matrices of the partonic processes qq — > ti r y and 
gg — > ti'j , and and M^ g are the NLO matrices for the q — q and g — g annihilation processes, 
separately. 

The virtual correction parts of the cross sections containing soft / collinear IR singularities. As 
we can see later that the soft /collinear IR singularities are canceled exactly after combining the 
virtual corrections to the pertonic processes pp(pp) — ► qq(gg) — > tij + X with the real gluon 
emission corrections and the gluon part of the PDF counterterms &G q 9 } g \/p (P — PiP j Q = 

16 




Figure 6: The pentagon Feynman diagrams for the partonic process gg — > tt'y . 
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u, u, d, d, s, s). 

D. NLO QCD corrected cross section for pp(pp) — > #7 + X process 

As shown in Eqs. fl2.6p . the PDF counterterms contain collinear IR singularities. By combining 
the contributions of the PDF counterterms with the hard collinear contributions of the qq —>■ ti'jg , 
gg — ► ti'-yg , q(q)g — > ti , ~fq(q) subprocesses, we get the expression for the remaining collinear 
contributions to the process pp{pp) — > ttj + X in 0(a s ) order as, 

a s T (1 — e) / 47T^ X ' 1 



do - 



coii 



27rr(l -2e) 



ij JO 



1 ,1 

dx\ / d%2 







x { [^i/Pi (xi, A*/) Gj/p 2 ( x 2, fJ-f) + Gi/Px Oi, A 4 /) <5i/p 2 (^2, A 4 /) 

a=i,j 
+ j) } 



(2.26) 



where 



AT(q -> g^) = C F (21n5 s + 3/2), -> ^) = 2iVln5 s + {11N - 2nj/)/6, (2.27) 



and 



G a / P (x,fif) = 



1— $a8„„i 



dy 



Ga'/p(x/y,fif)P aa '(y), 



Paaiy) = P aa >{y)\n ( 4 ) - KAvl 

y fjj 



(2.28) 



(2.29) 



where N = 3 and riif = 5, respectively. The explicit expressions for P aa i and P^ a , can be found 
in Ref . [24] . By adding the virtual correction, the soft real gluon emission corrections and the 
remaining collinear contributions shown in Eq. (12.26 p . the soft and collinear IR divergences are 
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vanished. The final result for the total QCD correction(Acr < 3 c ' D ) consists of a three-body term 
and a four-body term, i.e., Aa® CD = Aa^ + Ao^ 4 -*. The three-body term can be expressed as 

ss,cc,gg 

Aa^ = I da coLl + y 

^— ' _ 1 + ( 



ij=uu,dd 



i.l 



X 



Aa^ 



dxt / dx 2 \ (daYj + da s gij ) G i/Pl (x u n f ) G j/P2 (x 2 ,Hf) + (i j) , (2.30) 



G i/Pl (x u iif) G j/P2 (x 2 , nf) dfZ + (i <-> j) 



and the four-body term has the form as 



ss,gg 



1 



dx\ \ dx 2 

+ °ij JO JO 



ij=uu,dd 

+ / dx\ I dx 2 
Jo Jo 

u,d,s i i 

+ / dx\ I dx% 

Jo Jo 

q=u,d,s 



G c/Pl (xi, fi f ) G- c/P2 Oa, A 4 /) ^ g , C c + (c <-> c 



Gg/ft (Xl, /i/) Gg/p 2 (X 2 , flf) Vgg + (q^ g) 



(2.31} 



where &^ cE (s = x\x 2 s) and <^H(s = xia^s) are the cross sections for the partonic processes 
cc — > tt7g and zj — > ££7$ (ij = uu, dd, ss, gg) in the hard and hard noncollinear phase-space 
regions at the colliding energy s = x\x 2 s respectively, and (q = u,u,d,d,s,s) represent 

the cross sections in the noncollinear phase-space region for the partonic processes ug — ► ti^u, 
dg — ► i£~7<i, % — > ifyu, dp — > tt7<i, — > tt7s and sg — > tfys, respectively. 

Finally, the QCD corrected total cross section for the pp(pp) —>■ ti^y + X process is 

a QCD = a + Aa QCD = (J + A(J (3) + A(T (4)_ (3 35) 

In adopting Eqs. (12.301) . (12.311) . and (I2.32p for the numerical calculation, we use the CTEQ6M 
PDFs. 



III. Numerical results and discussion 

In this section we describe and discuss the numerical results for the LO and NLO QCD corrected 
physical observables for the processes pp(pp) — > ttj + X . We take one-loop and two-loop running 
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a s in the LO and NLO calculations, respectively[6j. The number of active flavors is Nf = 5, and 
the QCD parameters are Af ° = 165 MeV and Af^ = 226 MeV for the LO and NLO calculations, 
respectively. We set the factorization scale and the renormalization scale being equal(i.e., fi = 
[if — /i r ) and take fi = fi = m t by default unless otherwise stated. Throughout this paper, we 
take m t = 171.2 GeV and a{m z )~ l = 127.918[6J, and set quark masses m u = = m s = and 
m c = 1.3 GeV which are the same as the input parameters used in CTEQ PDFs[23j. The colliding 
energies in the proton-(anti)proton center-of-mass system are taken as = 14 TeV for the LHC 
and ^/i =1.96 TeV for the Tevatron Run II. 

To distinguish the photon from the jets requires the angle between the outgoing photon and 
jet constrained in the range of 1 j et > 8^ u j et (m the center-of-mass system of proton-(anti)proton). 
In our calculation we assume the produced (anti)top-quarks are always tagged and not effected 
by the phase-space cut. During our numerical calculation, we applied a number of checks to our 
calculations. 

1. We have compared the numerical results of the LO cross section for the process pp —* 
dd -> tij + X with y/s = 14 TeV by using FeynArts3.4/FormCalc5.4 [201 EI] packages and 
CompHEP-4.4p3 program[2H], and applying the Feynman and unitary gauges, separately. In 
the partonic luminosity integrations we adopt the CTEQ6L1 PDFs. The results are 125.1(1) fb 
(CompHEP, Feynman gauge), 125.0(1) fb (CompHEP, unitary gauge), 125.1(1) fb (FeynArts, 
Feynman gauge) and 125.1(1) fb (FeynArts, unitary gauge). It demonstrates all these results are 
in good agreement. 

2. We checked the UV cancelation both analytically and numerically. The IR finiteness is veri- 
fied numerically after combining all the contributions from the virtual corrections, renormalization 
constants, real gluon/light- ( ant i) quark emission partonic processes and the collinear counterterms 
of PDFs at the NLO. 

3. We used both the LoopTools2.2 package [2T] and our in-house library to check the correctness 
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of the deductions of the tensor integrals and IR-finite numerical calculation of the scalar integral 
functions at a few phase-space points. Both packages are developed based on the same expressions 
given in Refs. [291 1301 I3~T] . but in different codes. The numerical results are coincident within the 
calculation errors. 

4. The independence of the NLO correction on the soft cutoff S s and collinear cutoff S c , 
were proofed. Figures El^a) and [3(b) demonstrate that the total NLO QCD correction to the 
pp(pp) —>■ dd —>■ tij + X process at the LHC does not depend on the arbitrarily chosen value of 
the cutoff S s within the calculation errors, where we take the photon transverse momentum cut 

> 20 GeV, 0-yjet > 3° [in the proton-(anti)proton center-of-mass system] and 5 C = 2 x 10~ 6 in 
adopting the TCPSS method. In FigJT^a), the three-body correction[Acr^ 3 \ see Eq. (12.301) ]. four- 
body correction [Acr^, see Eq. (12.3 II) ]. and the total QCD correction (Aa® CD ) for the pp(pp) — > 
dd — > i?7 + X process, are depicted as the functions of the soft cutoff 5 S with 5 S running from 
2 x 10~ 5 to 2 x 10~ 3 . The amplified curve for the total QCD correction, Aa® CD , is presented 
in FigJ7|(b) together with calculation errors. Figures [Hl^a) and MJo) show the independence of the 
total NLO QCD correction to the pp(pp) — > dd — > ti r y + X process on the cutoff S c where we take 
S s = 2x 10~ 3 , 9 yJet > 3° and p^ ] > 20 GeV. In FigEKb) the amplified curve for Aa QCD for 
the pp{pp) —* dd — > tt'y + X process is depicted. The verification that the total QCD correction 
Aa® for the pp(pp) — > dd — > ti'-f + X process is independent of these two cutoffs, not only 
demonstrates the cancellation of soft /collinear IR divergencies in the total NLO QCD corrections 
to the processes pp(pp) — *> dd — » ti r y + X , but also provide an indirect check for the correctness 
of our calculations. In further numerical calculations, we fix p^ > 20 GeV, 6^ et = 3°, 5 S = lO" 3 
and S c = S s /50, if there is no other statement. 

In FigJHl^a) and Figs fTOT a) we present the dependence of the integrated LO and the NLO QCD 
corrected cross sections on the renormalization/factorization scale(/i) at the LHC and Tevatron 
RUN II, separately. There we assume fi = fi r = [if and define /i = m t . We can see that 
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Figure 7: (a) The dependence of the NLO QCD correction parts to the pp{pp) — > dd — > ttj + 
X process on the soft cutoff 5 S at the LHC. (b) The amplified curve for the total NLO QCD 
correction A<t® cd to the process pp{pp) — > cfd — >■ tij + X in Figj7(a), where it includes the 
calculation errors. 
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Figure 8: (a) The dependence of the NLO QCD correction parts to the pp{pp) — > dd — > ttj + 
X process on the collinear cutoff S c at the LHC (b) The amplified curve for the total QCD 
correction Aa^ CD to the process pp{pp) — * dd — > tt"f + X in FigE^a), where it includes the 
calculation errors. 



although the curves for the LO and NLO cross sections have visible variations when the energy 
scale /i runs from 0.1m t to 3mt, the curves for NLO become more stable in comparison with the 
corresponding curves for LO. It demonstrates that the NLO QCD corrections reduce obviously the 
dependence of the cross section on the introduced parameter \i in the plotted fi/no value range. The 
corresponding total K-factor [X = <jQ CD /a^o], the NLO QCD K-factor from the pp(pp) — > qq — > 
ttj + X processes [K q g = H — 9 ~"'g LC , ~ — ]> the K-factor from the pp(pp) — > gg — > ttj + X process 
[Kgg = 1 H — — ] and the K-factor from the pp(pp) — > gq(gq) — > ttjq(q) + X processes [X g g = 
x + 2^„, d ,^ +^ 9g - J j are plotted in Figl^b) for the LHC and FigfJOKb) for the Tevatron. 

FiglH^b) shows that at the LHC the integrated NLO QCD corrections always enhance the LO 
cross sections except in the range of 0.1 < /i//io < 0.21, and the NLO QCD corrections from 
the pp — > gg — > ti'j + X and pp — > qq — > tij + X(q = u, d, s, c) processes counteract the other 
contributions in the region of 0.1 < fi/fio < 0.46. Figj9]^b) shows when we take \ij \i§ = 1, the NLO 
QCD correction to the pp — > gg — *> tt'y + X process at the LHC is much larger than the correction 
to pp — > qq — > ttj + X processes, while Figure [TOTb) demonstrates the NLO QCD correction to 
the pp — > qq — > tt^ + X processes at the Tevatron is larger than that to pp — > gg — > #7 + X in 
the vicinity of /x//io — 1- in the further calculations, we fix /i = /i = m t . 

In Table [TJ we list the numerical results related to the data in FigsJH^a,b) for the LHC and 
Figs JTOT a.b) for the Tevatron at the position of /i = fiQ. They are the results for the integrated LO 
and NLO QCD corrected cross sections, the total K-factor [K = a ^ L ° ] of the process pp(pp) — > 
ti'y + X , the K-factor contributed by all the NLO QCD corrections to the pp(pp) — > qq — > 

££7 + X (g = u,d,s,c) processes [if g5 = 1 H 9= "'g LO 9g — ]> the K-factor contributed by the 

integrated NLO QCD correction to the pp{pp) — > gg — > ££7 + X process [X 93 = 1 H — ^ ] and 
the K-factor contributed by the corrections to the pp(pp) — > gq(gq) — > ti , ~fq(q) + X processes 
[Kgq = 1 + ^ Q ] at the LHC and Tevatron. 

The LO and NLO differential cross sections of the transverse momenta for the top quark and 
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Figure 9: (a) The dependence of the LO and NLO cross sections on the factoriza- 
tion/renormalization scale at the LHC. (b)The total NLO QCD K-factor for the process^ = 
cr® CD / <Jlo]-> the NLO QCD K-factor from the pp — > qq — > tt'y + X(q = u,d,s,c) processes 

[K qq = 1 + q = u % qq - ], the pp ^ gg ^ tt 7 + * Process ee 1 + Aajf D /<7 LO ] and the 

- V s _ ( Aa® CD + Aa^ D ) 

pp — > — > + X, (g = w, d, s) processes [Kg q = H — ^ 22 — ] versus the energy 

scale at the LHC. 2 ^ 
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Figure 10: (a) The dependence of the LO and NLO cross sections on the factoriza- 
tion/renormalization scale at the Tevatron. (b) The total NLO QCD K-factor for the process(fT = 
AcjQ /(Jlo), the NLO QCD K-factor from the pp — > qq — > tt'j + X(q = u,d,s,c) processes 

y^ s . c (Act®- 00 ) - 

[K q q = 1 + q= „ qq ], the pp -> ^ -> «7 + X process [i^ = 1 + Aof° D /a LO \ and the 

pp — > — > W7 + X, (g = w, d, s) processes [Kgq = H — ^ 22 — ] versus the energy 

scale at the Tevatron. 26 



Collider 


0"LO 


a QCD 


K qq 


K 99 


K 9q 


K 


LHC 


2.141(2) (pb) 


3.256(6) (pb) 


1.014 


1.272 


1.272 


1.521 


Tevatron 


43.79(4) (fb) 


42.80(6) (fb) 


1.148 


1.009 


1.009 


0.977 



Table 1: The LO and NLO cross sections, the total K-factor(.K" = ^ — ), the K- 
factor contributed by the NLO QCD correction to the pp{pp) — > qq — > ttj + X (q = 

u, d, s, c) processes [K q q = 1 H i=2h^ — n — j ; the K-factor contributed by the NLO 

Aa QCD 

QCD correction to the pp(pp) — > gg — > tt^y + X process [K gg = 1 H — ^ — ] and the 
K-factor contributed by the corrections to the pp(pp) — > gq(gq) — > ti^q^q) + X(q = 

u, d, s) processes [K gq = 1 + g p + g * j ] with /i = m t , # > 20 GeV and 

0«* et = 3° at the LHC and Tevatron. 



photon at the LHC, are depicted in Fig JTiT a) and FiglTTTb). respectively. The analogous plots at 
the Tevatron are depicted in Figs fT2T a) and (b). Figures [TIT a) and ITITb) demonstrate that the 
NLO QCD corrections enhance significantly the differential cross sections of pj) and p^! for the 
LHC, but the corrections make only a small impact on the distributions of pj) and p^r for the 
Tevatron as shown in FigsJTWa.b). In FigsJTWa.b) the NLO corrections at the Tevatron suppress 
the LO distribution of little bit except in the range of 20 GeV < pf < 100 GeV, while the 
corrections slightly reduce the LO differential cross section of p^ in the whole plotted p^! range. 
From the distributions of pj) in both FigJTTTb) and FigfT2Tb). we can conclude that most of the 
photons in the events of pp(pp) — > ti'y + X are produced in low transverse momentum range at 
the LHC and Tevatron. 

We adopt the definitions of the LO and NLO top-quark charge asymmetries in Ref.|32j, i.e., 

At _ a LO At _ °LO ( -i | ^ a NLO ^ a NLO \ (a i \ 

-^FB,LO — A FB,NLO ~ ~ ~ I 1 = + ] ' \° ' 

°LO a LO \ °LO a LO J 

In Eq.f l3.ip the notations a^ Q have the explicit definitions as 

Olo = °Lo(yt > 0) ± a LO (y t < 0), (3.2) 
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Figure 11: The LO and NLO distributions of the transverse momenta of the top quark and photon 
taking fi = m t , p^P > 20 GeV and 0^ et = 3° at the LHC. (a) for the top quark, (b) for the photon. 
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Figure 12: The LO and NLO distributions of the transverse momenta of the top quark and photon 
taking /j, — m t , > 20 GeV and 0^ et = 3° at the Tevatron. (a) for the top quark, (b) for the 
photon. 
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piUGeV) 


ctlo (fb) 


aQ CD (fb) 


A FB , LO (%) 


Afb,nlo{%) 


20 


43.79(4) 


42.80(6) 


-17.24(7) 


-11.41(8) 


30 


27.88(2) 


26.94(4) 


-18.21(8) 


-11.52(9) 


40 


19.09(1) 


18.26(2) 


-18.85(8) 


-11.90(9) 



Table 2: The LO and NLO cross sections and the LO and NLO forward-backward 
charge asymmetries of the top quark at the Tevatron with \x = m t , 0~ u f e t — 3° (in the 

proton-antiproton center-of mass system) and p F \ ut = 20 GeV, 30 GeV, 40 GeV, 
respectively. 

where cross sections (JLoiVt > 0) and (J LO {y t < 0) get the contributions from the top-quarks 
in the forward and backward hemispheres at LO, respectively, [The forward direction is defined 
as the orientation for incoming proton(Pi).], Aa^ LO denote the NLO QCD contributions to the 
cross sections cf . By using our program with the same conditions as used in Table 3.1 of 
Ref . |32j . we calculated the LO cross section and the top-quark charge asymmetry for the process 
pp -> ti + jet + X at the Tevatron. We obtained a LO = 1.582(2) pb and A FBLO = -7.70(6)% by 
taking PT,jet,cut = 20 GeV and fi = m t , which are coincident with those in Table 3.1 of Ref.[32j. 

In Table [2] we present the results of LO and NLO cross sections(o"Lo and <jQ CD ) and LO 
and NLO forward-backward charge asymmetries of the top quaxk(A t FB LO and A FB NLO ) for the 
process pp — > tij + X at the Tevatron. There we set fi = m t , 0^j et = 3°, and the photon 
transverse momentum cut p F l u t = 20 GeV, 30 GeV and 40 GeV respectively. The numerical 
results in the table show the NLO QCD correction reduces evidently the absolute values of LO 
asymmetry, e.g., in the case of p F \ u t — 20 GeV, the LO asymmetry |-A^ S LO | = 17.24% is cut 
down to \A FB NLO \ = 11.41% by the NLO QCD corrections. The absolute values for both LO and 
NLO asymmetries are quantitatively increased a little bit with the growing ofp^la from 20 GeV 
to 40 GeV. 
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IV. Summary 



In this paper we calculate the complete NLO QCD corrections to the top-pair production associ- 
ated with a photon at the LHC and Tevatron Run II. We investigate the dependence of the LO and 
NLO QCD corrected integrated cross sections on the factorization/renormalization energy scale. 
We present also the predictions for LO and NLO QCD corrected charge asymmetries of top-quarks 
at the Tevatron, and the LO and NLO differential cross sections at the LHC and Tevatron. We 
find from our numerical results that the NLO QCD radiative corrections obviously modify the LO 
charge asymmetry of top-quark, integrated and differential cross sections. And the uncertainty of 
the LO cross section due to the introduced unphysical energy scale /i, is significantly improved by 
including NLO QCD corrections. Our numerical results show that by taking /i — m t the K-factors 
of the NLO QCD corrections at the LHC and Tevatron RUN II are 1.524 and 0.977, respectively. 
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